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& ' Abstract 



Electrons in strong magnetic fields can be described by one-dimensional 
models in which the Coulomb potential and interactions are replaced by 
regularizations associated with the lowest Landau band. For a large class of 
models of these type, we show that the maximum number of electrons that 
can be bound is less than aZ+Z f(Z). The function f(Z) represents a small 
non-linear growth which reduces to A p Z (log Z) 2 when the magnetic field 
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B = 0(Z P ) grows polynomially with the nuclear charge Z. In contrast to 
earlier work, the models considered here include those arising from realistic 
cases in which the full trial wave function for N-electrons is the product 
of an iV-electron trial function in one-dimension and an antisymmetric 
product of states in the lowest Landau level. 



1 Introduction 

It is well-known that systems in strong magnetic fields behave like systems in 
one-dimension, i.e., a strong magnetic field confines the particles to Landau orbits 
orthogonal to the field, leaving only their behavior in the direction of the field 
subject to significant influence by a static potential. Motivated by this general 
principle and the work of Lieb, Solovej and Yngvason ^2] (LSY) on atoms in 
extremely strong magnetic fields, Brummelhuis and Ruskai jS] initiated a study of 
models of atoms in homogeneous strong magnetic fields in which the 3-dimensional 
wave-function has the form 

*(ri, r 2 . . . r n ) = ...x n ) T(y h z u y 2 , z 2 , . . . y n , z n ) (1) 

where T lies in the projection onto the lowest Landau band for an N-electron sys- 
tem. We follow the somewhat non-standard convention of choosing the magnetic 
field in the x-direction, i.e, B = (B, 0, 0) where B is a constant denoting the fields 
strength, in order to avoid notational confusion with the nuclear charge Z. 
The Hamiltonian for an N electron atom in a magnetic field B is 

N r 7i i 

H(N,Z,B) = J2 \ P > + A \ 2 -wTl TT ( 2 ) 



rj-r k \ 



where A is a vector potential such that V x A = B. The ground-state energy of 
H(N, Z, B) is given by 

E (N,Z,B)= inf (V,H(N,Z,B)V) (3) 
||*||=i 

Let EQ° nf (N,Z,B) denote the corresponding infimum restricted to linear combi- 
nations of functions of the form ((TJ. For extremely strong fields, it was shown in 
[12] that E /E™ ni -> 1 as B/Z^ 3 -> oo with N/Z fixed. 

In this paper we consider Ej(N,Z,B), the infimum when ® is further re- 
stricted to those functions of the form corresponding to a particular choice 
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for T. As discussed in is straightforward to show that 



Ej{N,Z,B) = v^B inf h(N, Z, B- 1/2 )<$>) + NB (4) 
ll*ll=i 

where 

- r 1 d 2 ~ i ~ 
h r (N, Z, M) = I ~ - ZV^x^l + Wj k (x 3 - x k ), (5) 

and we have scaled out the field strength B so that the only remnant of the mag- 
netic field is in the "mass" M = B- 1 ' 2 . The effective one-dimensional potentials 
Vj and Wj k can be written in terms of the functions HZj 

i r° s 2m e~ s2 



V m (x) = — r / -.sets, (6) 

which are discussed in Section EP1 and studied in detail in ^H]- The precise form 
of and Wj k depends on the choice of T; some special cases are discussed in 
Section El When T is a simple product of one-particle Landau states, or a finite 
linear combination of such products, the effective potentials satisfy V^ T (x) < V A 



x 



and Wj k {x) > V v {x) for some integers /x and v which depend upon T. We 
will primarily be interested in the case of symmetrized and anti-symmetrized 
products of one-particle Landua states. In this discussed in more detail in 

Section and Appendix El bounds of the form above are readily obtained. For 
those situations in which v < 2/x, a bound on the maximum negative ionization 
is given in Theorems and El 

In j3j we considered the simple, but unrealistic situation in which T is a prod- 
uct of Landau states with m = 0. In this paper we introduce a more realistic 
model, which we call the "Slater model", in which T is an antisymmetrized prod- 
uct of Landau states. As in j3] , we concentrate on the question of maximum nega- 
tive ionization. Define iV max (Z, B) as the maximum number of electrons for which 
the Hamiltonian (jSj) has a bound state in the sense Eq(N, Z, B) < Eq(N—1, Z, B). 

LSY ^2] showed that in extremely strong magnetic fields, atoms bind 2Z 
electrons in the sense 

Hminf ^M> 2 . (7) 

Z,B/Z3->oo Z 

and conjectured that 2Z was also an upper bound to this limit. However, even for 
the simple model in [3] we were only able to show the weaker bound N mauX (Z, B) < 
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2Z + l+cVB. Unfortunately, when B > 0(Z 3 ) as required for the limit in (J7J), the 
term c\fB = cZ 3 ! 2 dominates so that we can only conclude that N max (Z, B) < 
2Z + 0(Z 3 / 2 ). 

In this note we show that for a large class of one- dimensional models, including 
some in which T in (|1J) is a simple product or a Slater determinant, the bound on 
N mayL (Z,B) can be improved to one of the form aZ+Zf(Z) with /(Z) = O(logZ) 2 
when B grows polynomially with Z. Our results were announced earlier in [7j 
in the form 

N raax {Z,B)<aZ + Zg{Z,B) (8) 

with g(Z, B) = (log Z) 2 +log Z (log B) 1+UJ for some u> > 0. Subsequently, Seiringer 
[T%] gave a similar bound for the full 3-dimensional Hamiltonian. For fermions, 

Seiringer's bound [TH] has the form (jHJ) with g(Z, B) = min { (Js) 2 ^ 5 , ( log ]pr) 2 }; 
it is obtained by applying Lieb's method to the full 3-dimensional Hamiltonian. 
Hainzl and Seiringer [Oj then extended this bound to a density matrix model 
in which the variable perpendicular to the field is replaced by discrete angular 
momentum quantum numbers. 

Our use of one-dimensional models was motivated by a desire to understand 
the physics associated with the consequence of the one-dimensional character of 
atoms in strong magnetic fields, which is well-known and made precise in the 
work of LSY ^2]- Brummelhuis and Duclos [21 EH E] also showed that, for each 
fixed total angular momentum in the field direction, the full QM Hamiltonian (J2J) 
converges in norm-resolvent sense 1 to the projected Hamiltonian H S (N, Z, B) = 
U S H(N, Z, B)U S where U s denotes the orthogonal projection onto the lowest Lan- 
dau band. In the special case of zero total angular momentum, (^H S (N, Z, B)ty) 
has the form (j3J) when \l/ has the form (JTJ. Full details will be given in j^j. In 
contrast to [12] , the strategy in [21 E] does not require an a-priori bound on N in 
terms of Z, but does need to fix the total angular momentum in the direction of 
the magnetic field. 

Despite Seiringer's result ^Hj and the work in [Oj, we feel that our argument, 
which uses the RS localization approach, is of some interest. Because our analyses 
of the models in [H] showed that electrons are highly delocalized in the direction 
orthogonal to the field, it may seems surprising that such a localization technique 
works at all. However, a careful examination of the proof in section 3, shows 
that it reflects this derealization in the sense that the "inner ball" grows with 

: To be precise, let R be the resolvent of H, and R s that of H s . Then ||n s i?n s — R s \\ = 
0(B^ 1 /' 2 ) on T1 S H, and \\R\\ = OiB- 1 /' 2 ) on Tl^H at a distance (logS) 2 to the spectrum of H s . 
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B. Nevertheless, the localization error can be controlled with a modest excess 
charge as described in the results which follow. 

We now summarize our results in the theorems which follow. 

Theorem 1 Suppose the potentials in the Hamiltonian J3J) satisfy 

V J r (x)<V,(x) and Wj k {x) > ^K-r (-^) , (9) 

for all j, k and < v < 2 fx. Then for every a > there is a constant A a > 
and constants ai,a 2 (independent of a) such that the Hamiltonian h (N,Z,B) 
has no bound states provided that 

N >2Z + A a Z 1+a , and (10a) 
a ie za/4 > B > a 2 Z~<% (10b) 

where the exponent a can be arbitrarily small and the exponent 7^ depends upon 
v. In particular, when v = 0(1), it suffices to take 7^ > 2; when v = O(N) one 
must choose 7^ > 3. 

Although the non-linear term is higher order than 2Z, it is useful to write the 
linear term separately. It is due to the relative strength of the potentials near the 
nucleus, while the non-linear terms are needed to control the localization error. 
The upper bound on B is needed for technical reasons associated with the fact 
that the localization error can not be controlled when B grows exponentially. As 
discussed in Remark 1 of Section 13.41 the requirement v < 2\i can be relaxed at 
the expense of replacing 2Z by cZ in (jlUaj) with c > 2. 

The non-linear term in the lower bound ()10a|) can be improved to one that is 
logarithmic. We first state it in general and then under the simple, and realistic, 
assumption that B = Z p for p > 3. The case j u > 3 in (jllbjl . corresponds to the 
superstrong region considered by LSY in [T2] . 

Theorem 2 Assume that the potentials Vj and Wj k satisfy (0) with < v < 2\i. 
Then there are positive constants A, a e , a 2 such that the Hamiltonian h T (N, Z, B) 
has no bound states provided that 

N>ZZ + \ + AZ\ogZ\\og^\, and (11a) 
a e e zl/2 ~ e > B > a 2 Z^% (lib) 

where e > can be arbitrarily small and / y u is as in Theorem^ 
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Corollary 3 Assume that the potentials Vj and Wj k satisfy (GJ) with < v < 2\i 
and that B = aZ p for some a > and p > 3. Then there is a constant A such 
that the Hamiltonian h T (N, Z, B) has no bound states provided that 

N > 3Z + AZ (log Z) 2 . (12) 

Each of these results, immediately yields an upper bound on iV max (Z, B) which 
we state for ease of comparison with the results in [7113 CHI- Under the hypotheses 
of Theorem |21 

N max (Z, B) < 3Z + A \ogZ | log f | (13) 

Under the hypotheses of Corollary El 

N m UZ, B)<3Z + AZ(\og Zf. (14) 

Unfortunately, unlike Theorem^ the 3Z in the linear term includes a contribution 
from the localization error as well as the expected 2Z from electrostatics. 



2 Effective potentials 

2.1 Regularized Coulomb potentials 

The Landau state with energy B and angular momentum — m can be written 
compactly using the complex variable £ = y + iz as 

7£(l/, z) = [vrnr l,2 B^I 2 Z m e- B ^ /2 . (15) 

The effective one-dimensional potentials in our models can be written using the 
regularization of the 3-dimensional Coulomb potential with a Landau state, i.e., 

v£(x) = <7* n^>= / l -rr d y dz 

\ r \ JR 2 \ r \ 



roc g2rrig—Bs 2 



s ds 



1 

ml 



ml Jo Vx 2 + s 2 

00 ll m p -u 

du (16) 



o a/x 2 + u/B 



2B 



m+l 



ml 



e ^ 2 / {t 2 - x 2 )™ e - Bt 'dt, 

\x\ 



where r in R 3 and s = y 2 + z 2 . In view of the rescaling in it suffices to consider 
only the case B = 1 for which we drop the superscript, i.e, V m (x) = V^(x). The 
properties of V m (x) were studied in detail in [TH] and summarized in [7j. Those 
which we need here are listed below. 



V m (x) is monotonically decreasing for x > 0. 

1 



V m+1 (x) < V m (x) < 



\x\ 



(17) 
(18) 



> V m (x) > = 
V x 2 + m + 1 



\/ x 2 + m 
1 N ^ 

If Kv(x) = - J] V^x), then \4 v (x) < 2V N (x) 



j=0 



(19) 
(20) 



2.2 Simple product Landau model 

We restrict to wave functions of the form (JTJ) with T = Ylk=i ImAVk, z k) a simple 
product of Landau states. Then 



N 



k=l 



and 



where we rescale as in (j3J) and m = (mi, • ■ ■ , m n ). Then 

N 



h™(N,Z,M) = J2 



1 d 2 

Mdxl~ ZVm ^ Xj) 



+ ^W mjtmh (\xj-x k \), 
j<k 



V m is given by © (with B = 1) and the effective interaction satisfies 

Wm.m'Oc-z') = (Tm®7m')l ^7 7m ® 7m') 

r — r 



m+m' 



x — X 



3=0 
1 



v/2 J V ^2 
x — x/' 



(21) 



[V mi ... mjn H{N, Z, B)* mi ... mir ) = VB ($, h m (N, Z, B- 1 ' 2 )®) + NB (22) 



(23) 

(24) 
(25) 



for some bj > with ^2jbj = 1. That W mtm > can be written as a convex sum as 
in ()25|) was shown in |13| : For completeness, a proof in the special case m = m' 
is included in the Appendix. When m = m' = (J23|) reduces to W^oG^ — = 

72^°(^7r9 as snown in 

When all rrij = m are equal, we denote the effective Hamiltonian in (|2*H|) by 
h m (N, Z, M) and refer to it as the m-momentum Landau model. For this model, 

V r (x) = V m (x) and W r (x) = ^Y^my^j- The case m = was considered in 

and [7j. 

Corollary 4 Let h m (N, Z, B~ l l 2 ) be the Hamiltonian described above. Then for 
any a > and any B satisfying aie za > B > a 2 Z 2+£ for suitable constants a\ 
a 2 and some e > 0, there exists a constant A a > 

N max (Z,B)<2Z + A a Z 1+a . (27) 

Moreover, when B = aZ p for some p > 3, one can find a constant A (depending 
on a, p) such that 

N max (Z, B)<3Z + AZ(log Zf. (28) 



The Thomas-Fermi theories introduced by LSY in ^2] for the superstrong and 
hyperstrong regimes have a kinetic energy term typically associated to bosonic 
systems. It seems therefore reasonable to consider h m (N, Z, B~ l l 2 ) with domain 
in the symmetric wave functions. This is however in clear contradiction with the 
fact that the electrons described by the original 3-dimensional Hamiltonian (J2J) 
are fermions, and (JJ) should be anti-symmetric. Therefore, in the next section, 
we introduce a model which reflects the anti-symmetry. 



2.3 A Slater determinant Landau model 

It is reasonable to expect that the electrons will try to satisfy the Pauli principle 
by going into different orbits in the lowest Landau band, and that any realistic 
one dimensional model will have similar behavior. We now consider the special 
case in which \1/ has the form with T an anti-symmetrized product constructed 
using m = 0, 1, 2 . . . N—l. Thus, we let 

T= (— !=7oA---A7tf-i ) , (29) 
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where the wedge A denotes the anti-symmetric product so that T is a Slater 
determinant in the Landau states jj for j — 0, 1 . . . N — l. In this case, 



(H(N, Z, B)V, ^} = \ r B (h dct (N, Z, B- 1/2 )$, $) + NB 



with 



N r 

h det {N,Z,M) = 

3=1 



where 



1 d 2 

Mdxt~ ZVay{Xj) 



1 N ^ 
V„{x) = - V 3 

3=0 



+y^w det (\xj - x k \), 
j<k 



and the effective interaction is 



N-2 



(30) 



(31) 



(32) 



(33) 



3=0 



N-2 



with &2j+i > and b 2 j+i = 1. It then follows from (fTTj) that 



3=0 



(34) 



To verify (|32j) recall that the 7j's are normalized and mutually orthogonal. 
Thus 

1 / 1 \ 

— (7o A • • • A 7jv-i, -, — r 7o A • • • A 7^-1) 



N 



1 N ~ l 1 

k=0 ' 

-^14(x) = Kv(^) <2Vat(x) 



(35) 



N 



k=0 



where we used fl2TJ|) in the last line. The expression (jHB^) for W / det(a ; ) is proved in 
the Appendix. 
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Corollary 5 Let h det (N ', Z ', B 1 ^ 2 ) be the Hamiltonian described above. Then for 
any a > and any B satisfying a\e za/A > B > a 2 Z i+e for suitable constants a\ 
a2 and some e > 0, there exists a constant A a > 

N max (Z,B)<4Z + A a Z 1+a . (36) 

Moreover, when B = aZ p for some p > 3, one can find a constant A (depending 
on a, p) such that 

N max (Z, B)<6Z + AZ{\og Zf. (37) 
2.4 Other models 

For any fixed choice of Landau functions 7 TOl . . . 7 mjv the effective potentials V T 
and W r can be computed explicitly for both the case of a simple product and that 
of a Slater determinant. However, general formulas 2 are not so easily obtained. 
Nevertheless, the results in the appendix hold rather generally in the sense that 
W r is a convex combination of the form Yli=o ^775^ (yl) with J ~ 2max^ \mk\. 

Thus, one might hope to obtain bounds on the effective potentials similar to 
those in © but without the constraint v < 2/i. In fact, one easily finds W r (x) > 

^^-if^f) with v = 2max|mfc|. However, bounds better than V r (x) < V\{x) 
are not so easily obtained. In situations in which bounds of the form 

V r (x) < cV, (N) (x), W r (x) > ^V u{Nyi (-j=) (38) 

hold with the dependence of /i and v on iV known, this would lead to similar 
bounds on iV max (Z, B) with the contributon of 2Z to the linear term replaced by 
one of the form kZ with k depending on the relative size of n{N) and ^(iV). 

3 Proofs 

3.1 Localization 

The proof of the main theorem will use the RS localization method which is 
summarized in 0. The argument used here requires some refinements discussed 
in more detail in |T%1 ITB] . 

2 In the special case, T ~ JIaLi 7m fc (yfc' z fc) with all odd, W r is given by a convex sum 
which contains only V<y with even subscripts as in (|74ll . 
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Let Go, G\, ■ ■ ■ ,Gn denote a partition of unity consisting of functions which 
are Lipschitz continuous on and satisfy Ylj=o ^1( x ) = 1 as we U as the follow- 
ing additional properties: 

(i) supp(Go) C {x : \\x\loo < (1 + 5)p}, 

(ii) supp(G fc ) C {x : WxW^ > p, \x k \ > ^H^Hoo} for 1 < k < N, 
(hi) YlJ=o |VGj| 2 < A ^°f on supp(Go), and 

r s |2 ^ . (logiV) 2 (logAO^ 

(iv) E,=o I VG,| 2 < A §2 ^ < A on su PP( G fc), 

where A is a constant and V denotes the gradient in all the variables X\, . . . , xn- 
The existence of a partition with these properties is guaranteed by the construc- 
tions in [T^ITH]. 

In many applications, one wants a "sharp" localization which is achieved by 
choosing 5 so that 5 — > as Z — > oo, e.g, 5 = Z~ a for some a > 0. In such 
situation, it often suffices if ||x||oo < 2p on supp(Go)- In this paper, the term 
"localization" is a bit of a misnomer, as the radius p of the "inner ball" will grow 
with Z. In this case, it can be advantageous to let the localization be far from 
sharp and even permit 5 to grow with Z. 

Now let h(N, Z, B -1 / 2 ) be as in (jHJ, and note that the IMS localization formula 
[8 j implies that for any <&(xi, • • • , xn) in the domain of h = h(N, Z, B^ 1 ^ 2 ), 

N 

($,/>$> = ^2(G v $,hG u $)-($,LE(x)$) (39) 

t/=0 
JV 

= [h(N, Z, B- 1 ' 2 ) - LE(x)} G V <S>), (40) 

i/=0 

where LE(x) denotes the localization error 

N 

LE{x) = VB ^|VG,(x)| 2 . (41) 
It follows from properties (iii) and (iv) of G k that the localization error is bounded 
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VB(\ogN) 2 

L = A ^-^ on supp(Go) (42a) 

o z p z 



above by 



L k = \ v K ° / on supp(G fc ), k = l...N. (42b) 
d 2 p\x k \ 

To prove Theorem^ it suffices to show that 

(G fe $ [h(N, Z, B' 1 ' 2 ) - L k ] >e (N-l,Z,B) || G fc $ 2 1| (43) 

for k = 0,1... N with 6q(N, Z, B) the ground state energy of h(N, Z, B^ 1 ^ 2 ). 
Since (J4~3~)) is equivalent to 

[h(N, Z, B- 1 ' 2 ) -L k - e (N -1,Z, B)} > 0, (44) 

it suffices to show that the quantity in square brackets in (|4^1 is positive on 
supp(Gfc) for each k = 0, 1 . . . N. It is useful to handle the cases Go (inner ball 
estimates) and G\. . .Gjy (outer estimates) separately. 

In the next two sections, we use the convention that c and C denote constants 
in the sense that some constant exists for which the indicated bound holds. 



3.2 Inner ball estimates: 

On the inner ball (i.e., on supp(Go)) the 1-dimensional Hamiltonian (jSj) with 
effective potentials satisfying can be bounded by 

h(N,Z,B-^) > Ne (l,Z,VB) + W£>V^(^) (45) 

> - CN ^ os§ r + ^ 7m y w (4 6 ) 

where we used h(N, Z, M) > Ne (l, Z, M) + Y. j<k ^jk(\ x j ~ x k\), \xj - x k \ < 
2(1 + 5)p on supp(Go), property (fT§|) and 

d 2 „ _ ^ „Z 2 



~ zv ^ x) - eo(1, z ' ^ - ~ c 7m ( log ^ ■ (47) 



The lower bound in (J47j) above follows from the asymptotic formula in PQ for 
the ground state energy of the one-electron hamiltonian on the left in (|47|). Now 
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(ignoring the difference between N and N—l), the right side of (jjBj) will be positive 
if 

n , , 1/ 1 N 2 B , . 

(1 + 5) V + - < z2 . 4 - (48) 



2 G ^ 4 (logf) 



Since we can assume Z < N, the right side of (|48J) will be greater than z/ for 
sufficiently large N, Z if 

5 > ^Z^ (49) 

for some e > 0. For bosonic models we are be primarily interested in the case 
v = 0(1) for which ()49|) holds if B > Z 2+e ; for anti-symmetric models, v = O(N) 
so that we need B > Z 3+e . For now, we assume that (}4"9*j) holds, in which case 
the requirement (J48j) can be rewritten as 

rx9 , 1 N 2 B v 1 N 2 B , , 

(1 + 5) V < T 4 - - < ^ 4 . (50) 

Cz 4 (log^) 2 ^Z 4 (log^) 

Thus we can ensure that the right side of (J4T?|) is positive by choosing 

1 nVb , . 

' = c TT*i^?f (51) 

for some constant c independent of N, Z, B, 5. 

Then, since -e (N - 1, Z, B) > 0, the condition flSJ) will hold for G if 



v/5 V b; 2v/4(1 + (5)V + 2z/ <*V 



The first two terms in (|52|) behave like +CiV-^i ( log and the third like 

— A(^jp-) 2 — — - v Ag^ f2 g B • Comparing these expressions, we find that (using 
the assumption that N > Z) control of the localization error requires 

^ 6 1 < (logiVP(logf) 2 ^) (53) 
We consider two cases of small and large 5 separately. 
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a) When 5 = Z a , the left side of (J53J) behaves like -4-. In this case (J53)l 
holds with p given by ()51|) provided that N > Z is sufficiently large and 
-B < Ce zl/2 ' for some e > a. 

b) When 5 > 1, we can use the fact that ^jp < 2 to see that ()53|) holds for any 
5 > 1 if iV > Z and B < Ce zl/2 ' for some e > 0. Alternatively, we can 
eliminate the upper bound on B by letting N grow with B as well as Z. In 
particular, (jo^j) holds for any 5 > 1 if iV > Z\og 

3.3 Outer ball estimates: 

For any k such that 1 < k < N, we can write 

h r (N, Z, M) = hj£(N-l, Z,M)-^^- ZV k r (x k ) + WT k ( Xj - x k ) (54) 

with the understanding that h k (N — l,Z,M) is the Hamiltonian obtained by 
omitting terms in (J3J) involving x k , but with potentials defined by the iV-particle 
state T. Let Ej ,k (N — 1, Z, B) denote the corresponding ground state energy. 
Since hJ(N-l, Z, M) > EQ' k (N-l, Z, B) and -£5 > 0, it follows that 

(G fc $, [h r (N, Z, B- 1 ' 2 ) -L k - E^' k (N-l, Z, B)] G k <$>) (55) 
> (G fc $, [-ZV k r (x k ) + J^Wjkilxj-x^-L^G^). (56) 

For simplicity, we henceforth omit indices j, k on Eq,V,W and assume 3 that T 
is a symmetrized or anti-symmetrized product. 

Now for x G supp G k we have that \xj — x k \ < \ + < (2 + <5)|xfc|, so that 

1 T/ (\ x j- x k\\^ 1 T/ /(2 + 5)|x fc 



W^-x k l) > ^(JS^i)>^( 



v/2 1,-1 V V2 ) ~ V2 V y/2 
1 _ 1 1 



3 Since the bounds on 1/ and are independent of j this is not a significant restriction. In the 
case of distinguishable particles, one need only compare to Ej' k (N — 1, Z, B) for a particular 
k. When indistinguishable particles are associated with an asymmetric product of Landau 
states, the full wave function must be a linear combination of states of the form Q with terms 
associated with irreducible representations of S n to yield the appropriate total symmetry; this 
is a more complex situation than the simple product model considered here. 
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where we used (fTTjl and (fI5|l. 

It follows from the upper bound in (|19|) that 



Thus 



V T (x k ) < V^x k ) < ) (5* 



zv T (x k ) + £ w r (x k ) > —;4t= + N 9 l i f = (59) 



If one ignores S, the right side of (|59jl is approximately 

-Z AT 1 



(60) 



which will be positive if N > 2Z and v < 2\i. 

This explains the origin of the linear term in Theorem It remains to take 
into account the localization error (|42bjl . When 5-^0 one need only choose 
iV — 2Z large enough to control (|42bJ) . When larger choices of 5 are made, one 
pays the price of an increase of ZS in the electrostatic estimates, as shown in ()64|) 
below. 

Substituting in and using the estimate ()42b|) . one finds that is 
bounded below by 

, T{Xk) ^=G k $), (61) 



(l + |) 2 ^ + f 



where 



The expression (161 j) will be positive, ensuring that (|44jl holds for A; = 1 . . . N, if 

> for |x| > p. Since > p on supp(Gfc), we find that, with p given by 

(ED, 

Z Z Z , 7 2 X 4 



2x\~ 2N B N' 



logf) 4 <l (63) 
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when N > Z and ({49}) holds. Thus, we can conclude that 
2 T(» > N-1-2Z 



x z + fi 

> N - 1 - 2Z(l + |) - 2A 



(1 + f) 2 x 2 + | (1 + £) 2 (logAQ 2 (logf ) 2 Z 2 



5 2 iV 



(logiV) 2 (logf) 2 Z 2 (l + 5) 



5 2 



(64) 



where the second inequality used the assumption v < 2/i. 

Now, we can analyze the iV-dependence of the right side of (}6"4*j) by writing it 
in the form 



2T( X )>N-^f Q -R 



(65) 



where Q, R are positive and may depend on Z,B,5 but are independent of N. 
The expression on the right is increasing in N. Hence, for any fixed choice of 
Z, B, 5, if it is positive for some critical N = N c , then it will be positive for all 



N>N C . 



3.4 Completion of proofs 

To prove Theorem^ choose iV c = 2Z + 1 + aZ 1+2a . Then 



2T(x) > Z 



' 7 2a [(l + 2a) logZ] 2 (21ogZ-logi?) 2 " 

aZ ~ c Wz^ 6 



Thus if 5 = 0{Z a ), then for sufficiently large Z, 



2T(x) > Z 



aZ 2a - c(log Z)\2 log Z - log B) 2 - Z~ a 



(66) 



(67) 



for 



which can be made positive for sufficiently large Z as long as logi? < C^f, 

some suitable constant C. This will be true if B < aie z " /2 for some constant a\. 
Replacing a by a/2 yields Theorem ^ 

To improve the Z a growth of iV to one involving only logarithmic terms, as 
in Theorem El one can not let 5 — > 0. Moreover, the term ^"ffi in (|54*j) implies 
that one can not decrease the localization error by choosing 5 large. Therefore, 
we simply take 5 = 1, in which case (jMj) becomes 



2T(x) > N 



2Z - Z -c 



(logA0 2 (logf) 2 Z 



N 



(61 
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We now set N = N c in (JUS) and write N c = 2Z + 1 + Zf(Z, B). Then 

2T{x) > Z f(Z,B)-l-c K S B) { S I SM ' . (69) 

L }{%) J 

The inner ball estimates in section 13.21 require B = OiexpZ 1 ^). Under this 
assumption, Theorem ^implies that f(Z, B) grows more slowly than Z a for some 
suitable a, so that log f(Z, B) < \ogZ for sufficiently large Z. Hence we can 
find a constant A such that the right side of (|fi9|) is positive when f(Z, B) = 
A\ log^-| logZ. (Note that | log(Z 2 / B)\ will stay bounded away from 0, due to 
the lower bound B > Z 1 ", / y u > 2, which ensures that f(Z,B) > 1 for big Z.) 
This proves Theorem El 

Since our hypotheses do not permit B to grow exponentially with Z and the 
case of greatest interest is polynomial growth, e.g., B = Z 3+£ , it is useful to 
restate our results under the assumption that B = Z p for some p > 0. In that 
case, we can conclude that there is a constant A, depending on p, such that (J69j) is 
positive for f(Z) = A(logZ) 2 , which proves Theorem El This also gives a bound 
of iV max (Z, ZP) < 2Z + AZilogZ) 2 . 

Corollary 3 follows immediately from Theorem |21 and the discussion above. 

To prove Corollary it suffices to observe that the hypotheses of Theorems 
and 121 hold with \x = m, and v = 2m. Unless m depends upon N, we now have 
v = 0(1) so that flU holds if B > Z 2+£ . 

To prove Corollary El for the Slater model, note that it follows from (|2~Uj) and 
(jS3that V r {x) = U av (x) < 2V N {x), and from flHIJ) that W T {x) > V 2N -z{x). 
Thus, h det (N,Z,M) satisfies that hypotheses of Theorems ^ and 121 with /x = N, 
v = 2N — 2, but with Z is replaced by an effective charge of 2Z. This has the 
effect of doubling the coefficient in the linear term and modifying the constant in 
the non-linear term. Since v is O(N), ()49|) requires B > Z 3+e . 

Remark: With more technical effort some of the hypotheses can be relaxed 
and/or estimates improved as sketched below. 



I. The condition v < 2/i is needed only for the bound 



< 1 

■ 2 +fi — 

implicitly used in (|64p. However, this is actually used only in analyzing 
the outer region for which one can assume \xk\ > p is much larger than 
0. Hence, with a bit more effort, this condition can probably be dispensed 
with. At worst, allowing v > 2\i would only change the coefficient of the 
(lower order) linear term. 
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The upper bound B < Ce zl/2 E is somewhat artificial. It arises because we 
have chosen to state our results in a way that emphasizes the dependence 
of N max (Z, B) on Z. More rapid growth of B will increase the confinement 
of the electrons in two dimensions, but make them more delocalized in 
the direction orthogonal to the field. Hence, it is not surprising that the 
localization error will be harder to control if B grows exponentially with Z. 

In the case of Theorem |2l one can eliminate the need for this upper bound 
in controlling the localization error on the inner ball by using the fact that 
f > log 5- Thi s 

was noted in remark (b) after (|53|). However, if B grows 
exponentially with Z, then the estimate log f(Z, B) < \ogZ used after (|p^j) 
will no longer be valid. The upper bound can be eliminated by allowing N c 
to grow sufficiently with B. The choice f(Z, B) = (log Z) 2 + \og Z(log B) l+lJJ 
for some uj > will suffice. This proves the result stated as Theorem 3.1 in 
[7j and given after (jHJ) in the introduction. 

The linear term in the Slater model is doubled because we use the estimate 
Kiv(^) — 2Vn(X) which gives an effective charge of 2Z rather than Z. 
Although this bound is tight near x ~ 0, it is used only in the outer ball 
where \xu\ > p and one would expect V w (x) ~ Vn(x) ~ Vn-i(x). (Note 
that since v = 2(N — 1) it would suffice to have a bound with p = N—l.) 

In fact, using results in [31111 one can show that 

2r 2 /I \ 
V w (x) = 2V N (x)-—( K --V N ^ 1 (x)j (70) 

"iV s 

Vx 2 + N ' "V^ 
When | a; | > p, and B > Z 3 this becomes 

1 . r,f(^zy 



2V N (x)-^== + 0(-). (71) 



V„(x) w 2V N (x)--== + o( 

\/x z -4- N V 



Vx 2 + N V N 
< V N (x)+C^ (72) 

for sufficiently large Z and N > Z 1+a Thus, one could expect to show that 
(}36|) can be replaced by iV max (Z, B) < 2Z + A a Z l+a in Corollary El by using 
more refined estimates for V av (x). This would give the expected behavior 
of 2Z for the linear term when B grows more rapidly than Z 3 . 
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4 Discussion 



The first step of the RS method is to divide the system into a small "inner" 
ball in which binding is precluded because the electrons are confined to a small 
region, and an "outer" ball in which the localization error becomes negligible as 
Z — ► oo. For bosonic systems, one expects to be able to squeeze the electrons 
closer together, yielding a smaller cut-off p than for fermions. This feature is 
the only factor which precludes extending the proof of asymptotic neutrality in 
[TTj to bosonic atoms. This suggests that the localization error is not simply a 
technical artifact, but a reflection of a real physical effect. In the one-dimensional 
models considered here, the anti-symmetry required by the fermionic nature of 
electrons is achieved entirely within the Landau band. This results in a in a one- 
dimensional model that is bosonic, with the anti-symmetry reflected only in the 
effective potentials 

The one-dimensional confinement also delocalizes the electrons. This is re- 
flected by the the effective mass of M — B^ 1 ^ 2 in (JEJ which implies that in 
strong fields the electrons behave like extremely light particles. The uncertainty 
principle then implies that trial wave functions which localize the electrons can- 
not yield bound states. Thus, it may seem rather surprising that localization 
methods can be applied successfully. For atoms in strong magnetic fields, this 
terminology may be misleading because the cut-off radius p is not small. Instead 
p ~ NVBZ- 2 (\ogf-)- 2 which grows with B. Thus, localization methods can 
be used to obtain (non-optimal) upper bounds on iV max despite the fact that the 
electrons are highly delocalized and the size of the "inner" region becomes large 
as B — > oo. 

Lieb's method can be interpreted as a different type of localization in which 
Gk is essentially the inverse square root of the potential. In three dimensions, 
the resulting localization error can be completely controlled by the kinetic energy, 
eliminating the need for an additional inner/outer deloclization provided that the 
magnetic field goes to zero at infinity. However, as discussed in section 3 of [5 , 
control of the localization error is more complex in models resulting from the 
types of magnetic fields considered here. 

When Lieb's method was applied to a one-dimensional model in [5] control 
of the localization error led to a y/B growth in N. In [TSj Seiringer showed 
that better bounds can be obtained if one applies Lieb's method to the full 3- 
dimensional Hamiltonian, yielding results comparable to those obtained here. 

One might try to combine Lieb's method with the inner/outer localization 
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2 



used here, i.e., in Section EH] use Go as in (i) but in (ii) replace Gk by ,, Y 

V * \%k) 

for k > 1. In the case of the simple 0- model, the argument given in Section 4 of 
j3] can be used in the outer region with the term lim ^j^j replaced by ■ 



With p as in Section 13. 2\ this would yield a net bound of 

N max {Z,B)<2Z + AZ \log§\ (73) 

which is a very slight improvement. 4 To extend this to m ^ 0, would require 
additional work. However, for \x\ > p one should be able to show that [Kn(a;)] _1 ~ 
\/x 2 + m to obtain a similar bound. In the case of the Slater model, one would 
also need estimates of the type discussed in Remark 3. 
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Technology. It is a pleasure to acknowledge the hospitable environment at these 
institutions. 



A Appendix 

We begin by proving a special case of Proschel, et al's result |T3] that the effective 
interaction W mm i[x\ — x%) can be written as a convex combination of potentials 
of the form -j^Vj (^^p) with j < m + m'. In the special case m = m', only terms 
with even subscript occur in the convex combination. 

Lemma 6 The effective interaction W mm (xi — x%) defined in \25\i satisifies 

m 

W mm ( Xl -x 2 ) = E^^i(^). (74) 

j=0 

with b 2 j > and Yl^jLo — 1 • 

4 The change from a quadratic to a linear dependence on \ogZ is due to the fact that the 
LE arising from Gq does not require a (logiV) 2 in the numerator in (hi) and (iv), resulting in 
a net bound of the form ^ . 
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Proof: Substituting for •y m in (jzljl and writing out the resulting integral yields 

W mm { Xl -x 2 ) = / / d^d^ ===. 75 

™ ! Jc Jc v /(x 1 -x 2 ) 2 + |Ci-C2| 2 

We now make the complex change of variables to 

CT = -*=( Cl + Ca ) r = ^( Cl _ C2 ) (76) 

» = H = ^ICi + C 2 | * = M = ^sICi - Cal 

and let 6* be the angle between a and r. Then |Ci| 2 + IC2I 2 = + 1 2 and 

|d| 2 |C 2 | 2 = (s 2 + t 2 + 2stcos6)(s 2 + t 2 -2stcos9) 

= s 4 + t 4 -s 2 t 2 cos2#. (77) 

Substituting in (|75|) yields 

2 f 50 , 2 r 00 2 /" 27r (s 4 + t 4 -s 2 t 2 cos2#) m , 
WW ^i -x 2 ) = — e"* tdt / e" s sck / ^ v ==/-. 78 

Wo 7o Jo v/(xi - x 2 ) 2 + 2t 2 

Performing the integral over s and 6 yields 



00 

2 



Pit 1 

W mm {x 1 -x 2 )= I e^tdt !l ' (79) 

JO (Xi - X 2 ) 2 + 2t 2 

for some polynomial P(u) of degree 2m. Writing -P(w) = ■> anc ^ substi- 

tuting in (|79*|l immediately yields 

2m 

W mm (xx-x 2 ) = ^-L^(^). (80) 

i=0 

It remains only to show that the coefficients \ are even, positive and sum to one. 
Applying the binomial expansion to the numerator in (|78|). one easily sees that 
terms with i = 2j have positive coefficients b 2 j > 0. When i = 2j + 1 is odd, one 
has an integral of the form J Q n cos 2 - ?+1 29 d6 = which implies b 2 j + \ = 0. Thus, 
only the coefficients with % = 2j survive, and these are positive. To see that 
Y^=ob 2 j = 1, it suffices to use the fact that both W mm (x) and all Vi(x) behave 
like l/\x\ at infinity. 
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When m m' the integrand will contain an additional factor of the form 
(s 2 + t 2 + 2st cos #)l m - m I. As above, one can show that terms involving cos- 5 9 
integrate to zero when j is odd. 

One can use a similar strategy to show that any inner product of the form 

(7i(Ci)7*(C), |^7;(Ci)7*(C2)} or (81a) 

(7i(Ci)7 fc (C 2 ), F 4^7.(Ci)7i(C2)> (81b) 

j+k 

can be written as a linear combination ^^Cj 4g Vj ( ( Xl ~^ ) with = 1. In the 

i= ° 

case of (|81aj) . as sketched above and shown in t ne coefficients c, are positive 
giving a convex combination. For the exchange integrals (j81b|) . this need not 
be true. However, for antisymmetric products, the two types of integrals arise 
together in combinations whose net coefficients are positive, as shown below. 

Lemma 7 For any choice of < ni\ < m 2 . . . < m^, 

W™t- mN (xi - x 2 ) = ^(7 mi A---A7 mN ,: r7mi A ■ ■ ■ A7 mjv ) 



i"! - r 2 



= ^EWV(^). (82) 

j=0 

with b 2 j+i > and ^ . b 2 j + i = I, and J = m N -i + m N , 

In the Slater model, we have = k—1, from which it follows that J = 2N—3. 



Corollary 8 For the Slater model, the effective interaction in h det (N,Z,B 1 ^ 2 ) 
J3J) is a convex combination ofV m with odd m — 1,3, .. .2N — 3, i.e., 



N-2 



W r (x 1 -x 2 ) = W dct ( Xl -x 2 ) = ^^b 2j+1 V 2j+ i(^j§*-). ( 83 ) 



3=0 

N-2 



with b 2 j + i > and b 2 j+i = 1 • 

3=0 
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Proof: To prove Lemma UJ we first consider the special case N = 2. Let j, k be 
fixed, and write 



— ^2) — *\ ij/y ik 



W^f (X! - x 2 ) = 2( 7i A 7fc (Ci, Ca), | r 7iA 7fc (Ci, C2) ) 

l r i _ r 2| 




c Jc 



V(^i-^) 2 + ICi-C 2 



rfCi^C 2 , (84) 



where we used (fH)]) and as before, r,- = (xj, yj, Zj) and Q = yj + izj. Now make 
a change of variables as in ([76)1 . Using the binomial expansion, we find 

C i c * = 2 -0+*)/a £(-1)"® Q r « + V+ fe -^\ 

i/=0 ^=0 

so that 

a=0 

where 

Aq = 2 -o-+*)/2 ^ [(-ir-(-ir])(9©- 

+ = a 
f < j, // < fc 

and we suppress the dependence of v4 a on j, k. When a is even, (— 1) M = (— l) v 
so that A a = 0. Therefore, only terms with a odd survive in the sum (|85|). 
Moreover, 



c ./c 



lC/C 2 fc -Ci^| 2 e-^' 2 -^' 2 

> /(x 1 -X 2 ) 2 + |Cl-C2| 2 



^(x!-x 2 ) = / / ^ SlS f„ '' =^0^ (86) 



a odd /3 odd 

Next, write r = ie*^ and use the fact that 



> > A a A / / . e 8 1 dadr 

"JcJc ^( Xl - X2 y + 2ti 



r*oo p2ty 

T a T P f{t)dr= I t a+p f{t)dt / e^-fidtp 
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is zero if a 7^ j3 to see that the integral after (jSfij) becomes 



\2a+l 



(2tt) 2 Y \A a \ 2 / s^+^+V* 2 ^ / ' ' = rffc. (87) 

Integrating over s then yields 

a odd 

for suitable constants b a , of which we only want to note that they are strictly 
positive when a is odd and in the range 1 < a < j + k. As before, bj = 1 
follows easily from the fact that both Wdet(x) and the Vj(x) behave like l/\x\ at 
infinity. This proves Lemma [3 in the case N = 2. 

The general case then follows from the fact that an iV-particle Slater deter- 
minant is a convex combination of two-particle slater determinants. In fact, 

W£b"' m "(xi ~ x 2 ) = ^(j mi A • • • A 7miV , 1 7mi A • • • A 7 mjv ) 

ri — r 2 



N(N 

j<k 



2 IjT 'y^ilrrij ^lm h (Cl) C2) ) , " f Imj Mm k (G, C2) ) ■ 

ir 1 — r 2 | 



2 



JV(jV-l) 

j<k 
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